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Abstract 

We start from the two existing QCD evolution equations for struc- 
ture functions, the BFKL and DGLAP equations, and discuss the 
theoretical hints for a unifying picture of the evolution in x and Q 2 . 
The main difficulty is due to the property of angular ordering of the 
gluon radiation driving the evolution and the cancellation of the re- 
lated collinear singularities. At the leading log 1/x and leading log Q 2 
accuracy, we find a unified set of equations satisfying the constraints. 

1 The two QCD evolution equations and the 
unification problem 

There is a well-known method to obtain the QCD-theoretical predictions for 
quark and gluon structure functions measured in Deep-Inelastic Scattering: 
the resummation of leading logarithmic contributions at all orders of the per- 
turbative theory; indeed, the existence of collinear and infrared singularities 
in the evaluation of radiation corrections to the point like lepton-parton scat- 
tering leads to effective coupling constants of order unity, and thus to the 
need of resummation techniques. Two types of resummation do exist. 
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At large Q 2 = —q 2 , where q is the quadri- momentum transferred to the 
target-proton, the collinear singularities lead to effective coupling constants 
of order a s (Q 2 ) In (Q 2 /Aq CD j . The resummation to all orders of the leading 
logarithms (LLQ 2 ) leads to the well-known "Altarelli-Parisi" (or DGLAP) 
evolution equations 1,2 - 1 . 

At small values of the Bjorken variable x = Q 2 /2p.q (p is the quadri- 
momentum of the target), a similar problem appears with the soft part of 
the gluon radiation, namely an effective coupling constant of order 0:5 In 1/x. 
It is thus necessary to resum the corresponding leading logarithms (LLl/x) 
to all orders. This non-straightforward resummation has been first performed 
by L. Lipatov and collaborators (BFKL) and leads to a singular behaviour 3 ) 
of structure functions at small-x. This result has been recently revived by 
Hera results on the quark structure function in the proton at very small- 
x. They have revealed a behaviour in qualitative agreement with this QCD 
prediction 4 ^ . 

For both phenomenological and theoretical reasons, it is interesting to 
address the problem of unifying the two mentionned equations into a single 
scheme. Since Hera experiments 5 ) cover a very large range in x and Q 2 , 
it is quite important to have a unified description of the QCD-evolution 
of structure functions in the whole a;— range. Moreover, it could solve the 
dependence on initial conditions for the evolution equations, which has to do 
with the unknown non-perturbative regime of QCD. 

On a more theoretical ground which is of concern in the present paper, 
it is to be remarked that the unification problem has already been suggested 
and discussed in the past. It has first been noticed that both LLl/x and 
LLQ 2 can be formally taken into account by a suitable combination 5 ) of 
the evolution kernels. On a more rigorous basis, it has been shown that a 
uniform description of the gluon radiation responsible for the evolution of 
structure functions in the whole x-range is possible due to the property of 
angular ordering 6 \ Within that picture, it can be shown that the collinear 
singularities present in all gluon production amplitudes contribute only in 
regions satisfying the following kinematical property: 

Q/x > ... > 6i > 0i_i > ... > 0i (1) 

where 9i xs (qr)i/xi are the angles of the emitted gluon with respect to the 
direction of the first emitted gluon momentum. One can separate two cases: 
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i) at fixed x, corresponding to finite (non strongly ordered) Xi, one recovers 
the well-known g r ordering of the LLQ 2 resummation technique 2 ^ . 

ii) at small x, the gluon momentum fractions Xi are necessarily strongly or- 
dered ( <1] and thus q r ordering is not implied by the relations (1) and 
angular ordering is expected to contribute to LLl/x singularities, together 
with the infrared singularities. 

However, it remains an important constraint to be fulfilled by any unifi- 
cation scheme based on angular ordering. The Lipatov equation is only re- 
covered provided that these LLl/x singularities related to angular ordering 
exactly cancel in the evolution of structure functions 6 ' 7 \ This cancellation 
is not valid for other observables such as multiplicities, average transverse 
momentum etc... As we shall see now, this stringent constraint leads to 
non-trivial consequences on the unified evolution equations. 

Let us first write the Altarelli-Parisi equations in a suitable form for unifi- 
cation. For this sake, we will restrict ourselves to the case of a fixed coupling 
constant a$ (as for the BFKL derivation) and consider the double inverse 
Mellin transform of the singlet (F s ) and gluon (F G ) structure functions : 

F s ,g^Q 2 ) = J^ er^'J ^ eWWtpsM-y)- (2) 

The Altarelli-Parisi equations (for fixed dig) can be written in matrix form 
for (ps and ipc as follows. 

as_ ( v F 2n F (j) F 
4tt 7 ^ 0£ 

where |z/ G , v F , 0g, are the usual (j-dependent) Altarelli-Parisi weights 2 ^, 

and (p^G are the initial conditions. 

Now, one has to modify the equation (3) to take into account the BFKL 
contribution in the gluon sector. The dominant contribution can be expressed 
3 ' 5 ) as a singularity in the j-plane situated at the value 

a s N c 

Jl = 1 + X (7) (3) 

71 

where 

X (7) = 2^(1) - 1K7) - ^(1 - 7); 1K7) = ' 1 l "' ) 
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is the eigenvalue-function of the BFKL kernel. Indeed, assuming a simple 
pole singularity tps,G (j — Jl) _1 and inserting it in the Mellin transform 
(2), one gets: 

FsgM) = f P~ e 7lnQ2/A2 e* S ^^l* ~ (^V**^ 2 , (4) 
J 2m ^ A ' 

where one has made use of a saddle point method to integrate around the 
point 7 C = 1/2, x(lc) = 4 In 2. Equation (5) corresponds exactly to the 
leading BFKL behaviour 3 ), up to logarithmic corrections (which would be 
determined by the nature of the singularity in the j-plane). 

A consistent modification of equation (3) in order to complement the 
singular behaviour (4) is the following 8 -* ; let us replace the gluon contribution 
to the anomalous dimensions: 

v G{i) — ► v*G(i) = 7 x{l) {vg + ~ *, (5) 

where ^ is an arbitrary function holomorphic in the j'-plane near j = 1 
and below. Such a modification inserted in equation (3) provides a formal 
unification of the Altarelli-Parisi and Lipatov kernels. Indeed, inverting the 
relation (3), after the replacement vq — > i/q, one gets 



with 




Now, the solutions of equation (6) depend on the region in the complex j- 
plane involved in the Mellin transform (2), and thus on the region in lnl/rr 
one is investigating; two cases appear: 

i) when x is not small, aslnl/rr 1, the modification (5) has no effect, 
since the zeroes of D(j, 7) are obtained for small values of 7 (of order as)- 
In that limit, one has from the very definition of x{l) '■ 

xirt) » 1/7 + 0(7 2 ); ~ vq + o{o?). (7) 
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One recovers the ordinary Altarelli-Parisi scheme 3 -* and the corresponding 
evolution equations (at fixed as)- 

ii) When a^lnl/x = 0(1), the singular structure of the BFKL kernel plays 
a role, driving the relevant domain of the Mellin integration over 7 near to 
the "critical" value 7 C = 1/2. 

In those conditions one recovers the singular behaviour compatible with the 
BFKL calculations. Taking the appropriate limit j — > 1, a~s/{j — 1) = : 

n( ■ \ 1 ® N cX(l) , a 4 In 2 

D{J, 7) a 1 - -r-r. — 7T. ~ 1--N C - (8) 

47T(j — 1) 3 — ►l 71 J — 1 

At first sight, the DGLAP and BFKL evolution equations can be unified 
for an arbitrary regular function ^f(j) in eq.(6). For instance let us consider 
the combination: 

v* G + is F = 7x(t) {»g + ^} + vf - * (9) 

appearing in -D(j, 7) at first order in as- Following the arguments of refs. 6,7 ^ 
as we have stressed upon in our introductory discussion, the quark-loop con- 
tribution up, which is present at fixed value of X clS db result of collinear 
singularities, should be absent from the evolution equations for small value 
of x. More precisely, if not cancelled appropriately, it would bring a new 
LLl/x singularity, due to the angular-ordering property including emitted 
quarks. It is thus compelling to choose when j — > 1 in order to 

obtain the desired cancellation. This is just the mechanism proposed in our 
paper 8 ). Indeed, the problem has been noticed to arise when one is to include 
"finite parts" into the evolution equations at small-a; 9 ). 

As a consequence, considering the proposed cancellation to be valid in 
the j-plane around the leading singularity jz,, one writes 

» Mjl) (10) 



aN c log 2 



n -1 



7T 



As noticed in Ref. [8], the equations (10) lead to an appreciable modification 
of the location j L of the BFKL singularity endpoint, in better agreement 
with phenomenological determinations 4 -*. 

Among other interesting properties, the set of equations (10) ensures (at first 
order in as) that the quark loops do not contribute to the small-x evolution. 
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Moreover, the expression of D(j, 7 ) preserves the position of the saddle 
point in the 7 plane at the critical value 1/2, as expected from the conformal 
properties of the BFKL kernel 10 ^. 

As a conclusion, the unification of the evolution equations for structure 
functions appears possible, at least in the leading logarithmic approximation 
and at fixed coupling constant as- Despite stringent constraints due to the 
mismatch of LLl/x and LLQ 2 perturbative resummations, a unified set of 
equations for the whole range in x and Q 2 can be written and leads to non- 
trivial predictions. A number of interesting questions remain open for future 
investigation, let us list some of them: 

1) Is it possible to implement unified evolution equations at the next-leading- 
order ? 

2) In the same context, how the result might be influenced by the running 
of as ? 

3) What are the phenomenological consequences of unified equations ? 

We hope to be able to provide answers to these questions in the near future. 
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